Storage of classical information in quantum spins 
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Digital magnetic recording is based on the storage of a bit of information in the orientation of a 
magnetic system with two stable ground states. Here we address two fundamental problems that 
arise when this is done on a quantized spin: quantum spin tunneling and back-action of the readout 
process. We show that fundamental differences exist between integer and semi-integer spins when 
it comes to both, read and record classical information in a quantized spin. Our findings imply 
fundamental limits to the miniaturization of magnetic bits and are relevant to recent experiments 
where spin polarized scanning tunneling microscope reads and records a classical bit in the spin 
orientation of a single magnetic atom. 
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Recent experimental breakthroughs have laid the foun- 
dations for atomic-scale data storage, showing the capa- 
bility to read and manipulate the spin of a single mag- 
netic atom with a spin polarized scanning tunneling mi- 
croscope (SP-STM)^, Read out is based on tun- 
neling magnetoresistance at the atomic scale [1]: for a 
fixed current spin polarization in the tip, the resistance 
is higher when the magnetic adatom spin is antiparallel to 
it. Spin manipulation is based on spin-transfer torque at 
the atomic scale [H, |3|: angular momentum is transferred 
from the spin-polarized tunneling electrons to the mag- 
netic atom. We are specially interested in cases where 
the magnetic atom is weakly coupled to the conducting 
substrate, which can be done thanks to a monoatomic 
layer of insulating material [4]. As a result, the spin of 
the magnetic atom is quantized, and can be described by 
a single spin Hamiltonian, as revealed by inelastic elec- 
tron tunneling spectroscopy (IETS)[2, 0, Hf, identical to 
that of single molecule magnets [6]. 

In this letter we address two fundamental questions 
that arise when considering magnetic recording in the 
quantum limit, i.e., the storage and readout of a classi- 
cal bit of information on a quantum spin. First, what 
is the role played by spin-parity in the readout and con- 
trol operations of a quantized spin?. The two physically 
different ground states required to encode the two logi- 
cal states of a bit appear only in the case of semi-integer 
spinQ for which quantum spin tunneling [gI, Q is forbid- 
den. We also show that current induced single atom spin 
switching is only possible for semi-integer spin. The sec- 
ond question is how can the magnetoresistive single spin 
readout be performed without disturbing the spin state? 
Here we study the problem of the back-action, akin to 
the quantum non-demolition [qI problem on a decohered 
qubit. 

The physical system of interest consist on a magnetic 
atom with quantized spin 5'[2', 0, Ej. The magnetic atom 
is probed and controlled by a SP-STM. The quantized 
spin of an atomic scale nano-magnet on a surface can be 



described with a single ion Hamiltonian [2I, 

^spin = DS^, + E{Sl - SD + gfiB^.B, 



(1) 



where D and E define the uniaxial and in-plane magnetic 
anisotropy. The eigenvalues and eigenfunctions of ([T]) 
are denoted by Em and |M), respectively. The above 
Hamiltonian accounts for the measured lETS in = 1 Fe 
PhthalocvaninefToj. S = 3/2 Cobalt adatomsjEI, 5' = 2 
Fe adatoms and S = 5/2 Mn adatoms[4]. 

The Hamiltonian of the total system features a sin- 
gle ion Hamiltonian exchange- coupled to the transport 
electrons 0, [HI, [III , = + '^s + Hspin + V, where 
Ht + 'Ws = J2x a ^cr(^)^A a^A,cr dcscribcs the tip and sur- 
face electrodes, with quantum numbers A = 7^) and cr, 
the momentum, electrode {r] = T, S) and spin projection 
(J along the tip polarization axis. We assume a spin po- 
larized tip with polarization Vt and a spin-unpolarized 
substrate S. The V term introduces interactions between 
tip, surface and the magnetic atom: 
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with a the Pauli matrices vector and S the magnetic 
atom spin. Equation (j2j) describes both spin-independent 
tunneling, described by the 7^^^^ term, as well as spin 
dependent processes, described by the 7a, a' term, where 
carriers can either remain in the same electrode, provid- 
ing the most efficient atomic-spin relaxation channel, or 
switch sides, which gives rise to spin-dependent tunnel- 
ing current. Neglecting the momentum dependence and 
considering the initial and final electrode of the scatter- 
ing events, there are 6 non-equivalent exchange integrals. 
We describe them as 7^^^, = Vr^Vr^'T^ and 7a, a' = TjVr^Vr^' ^ 
with the spinless Tq and spinfull Tj tunneling matrix el- 
ements, and two dimensionless scaling parameters that 
describe the strength of the tip-atom and surface-atom 
single-particle hoppings, vt and vs- 
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The effect of V on the energy levels is considered weak 
in the sense that it can be described within lowest order 
Fermi golden rule. We assume that the correlation time 
of the reservoirs formed by the electron gases at the tip 
and surface is short enough so that non-markovian effects 
are negligible [13]. The dissipative dynamics of the atomic 
spin described by Hspin^ under the influence of the dis- 
sipative coupling to the tip and substrate, is described 
in terms of a Bloch-Redfield (BR) master equation in 
which the coupling to the reservoirs is included up to sec- 
ond order in the coupling V: dtp = —^[HspimP] + 
with C the Liouvillian that accounts for the Kondo cou- 
pling V[13]. This equation describes the evolution of the 
diagonal terms in the density matrix, the occupations 
Pm = Pm,m^ as well as the off-diagonal terms or coher- 
ences, pM,M'- In steady state, the density matrix pM,M' 
described by the BR master equation does not contain 
coherences, and only the diagonal terms Pm survive [i3|. 

The relevant scattering rates can written in terms of 
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where e is some energy scale relevant for the process in 
question, a can be or J, and is the density of states at 
the Fermi energy in electrode rj. The elastic conductance 
has a contribution coming from the spinless tunneling, 

= 2e^ ^^^e^^^ ^ which plays no role in the remainder 
of the manuscript (e is the (negative) electron charge). 
From the experimental linear conductance we get that 
-ffs(lmeV) = I/e - 0.1 - 5.GHZ0. 

The spin-readout is based on a second contribution 
to the elastic conductance coming from elastic exchange 
between transport electrons and the spin of the mag- 
netic atom that gives rise to spin- valve term in the total 
conductance [3|]: 



(4) 



where {S) is the expectation value of the electronic spin: 



(5)=^Pm(^)(M|5|M). 

M 



(5) 



Thus, for finite tip polarization, the conductance is sensi- 
tive to the expected value of the atomic magnetic moment 
along z. Thus, if the quantum spin can be in two differ- 
ent spin states at zero applied field, ideally with {M\S\M) 
parallel and antiparallel to the tip moment, then a mag- 
netoresistive readout of a classical bit of information on 
a quantum spin is possible. 

We now discuss the necessary conditions for the exis- 
tence two ground states. First, D should be negative. To 
see this, we consider first the idealized case of a quantum 
spin with £^ = 0. The energy levels are Em = DS^^ 
with Sz = ±5*, ±(6* — 1) ... If is positive, the ground 
states doublet would have Sz = ±1/2 for semi-integer 



spin, which can give rise to Kondo effect [5[, or S'^ = for 
integer spin. In both cases the magnetic moment is zero. 
Second, the spin should be integer. Kramer's theorem [l^ 
states that, at zero field and with E ^ 0^ integer spin 
systems have non-degenerate spectrum, but semi-integer 
spins have, at least, a twofold degeneracy. These zero- 
field splittings can be interpreted in terms of quantum 
spin tunneling, which is suppressed for semi-integer 
Thus, the E term splits all the doublets of the E = 
spectrum only for integer S. Zero field splitting for in- 
teger spins has a very important consequence, which de- 
rives from the following general result. For zero applied 
magnetic field, the matrix elements {M\S\M) are zero 
for every non-degenerate eigenstate of Hspini'^^]- Thus, 
from Eq. ([5]) we get that, at zero applied field, it is 
impossible to have a net magnetic moment for integer 
spins. In contrast, for semi-integer an arbitrary small 
magnetic field along an arbitrary direction Cl will choose 
between the two ground states ^+ and resulting in 
{g±\S ' ^\g±) ^ 0. These two states provide the physical 
realization of the two logical states of the classical bit. 

lETS confirms this scenario for Fe {S = 2) and Mn 
{S = 5/2) on CU2N In both cases, D is negative 
{Dfq = —1.55 meV, Dmu = — 39/ieV). However, in the 
case of Fe, there is a single ground state due to quan- 
tum spin tunneling induced by E^ with a null average 
magnetization. In contrast, for Mn, with 5 = 5/2 the 
in-plane anisotropy does not lift the degeneracy of the 
ground state doublet for the Mn, see Fig. [TjL),b) 

The storage of information in D < semi-integer spins 
is limited by spin relaxation, originated both by elastic 
and inelastic processes. The later, addressed below, are 
exponentially suppressed when both bias and tempera- 
ture are smaller than the excitation energy. In contrast, 
the rate of elastic scattering between the two ground 
states g± due to coupling to the substrate, for semi- 
integer S is given by 



(6) 



a=x,y,z 



with kB the Boltzmann constant. The wave functions 
\g±) satisfy 
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2n 



(7) 



where n = 1,2,...,S' — ^, and Cn are dimensionless num- 
bers of order 1. We see that in-plane anisotropy enables 
the exchange assisted elastic spin flip^] \{g-\Sa\g-\-)\'^ cx 
[E / DY^~^ . Thus, elastic population scattering is sup- 
pressed as either S or D/E increase. Numerical calcula- 
tion (see Fig. [2^) yield lifetimes in the range of microsec- 
onds for S = 5/2 and \D\ = bE at 0.4K. 

Importantly, the coupling of the atomic spin to the 
conduction electrons kills the coherence between the two 



3 



ground states g±^ which satisfies the equation dtpg^^g_ = 
—^g+,g-Pg+,g-' The decohereuce rate ^g+,g_ contains 
contributions from both the non adiabatic terms, that 
imply population scattering, like in Eq. ([6|), and adia- 
batic terms, which do not involve energy exchange with 
the reservoir [13j. The substrate mediated rate for the 
adiabatic term reads as: 



9+,9- 



{9+\Sz\9+) 



\Sz\9-) 



(8) 



Thus, coupling to the electronic environment kills quan- 
tum coherence more efficiently in states with opposite 
magnetic moments, acting as a which path detector [l7| 
and favoring the magnetoresistive readout. To leading or- 
der in E/\D\, we have Tf^^g_ ^ 1ss{^bT)S'^. In contrast 
with population scattering, decoherence rate increases for 
larger S and is independent of E/\D\. Thus, larger S fa- 
vors spin memory but kills quantum effects. The ratio 
of the adiabatic decoherence rate, Eq. (|8]), and the elas- 
tic population scattering, Eq. ([6]), reads S'^ {\D\/ E)^^~^ ^ 
which is above 10^ for Mn in CU2N. 

We now turn our attention to the effect of parity on the 
process of magnetic recording, based on atomic scale spin 
transfer torque, which has only been studied for semi- 
integer spins so far 1, Current flowing through the 
spin-polarized tip, transfers angular momentum to the 
atomic spin. When the transfer rate exceeds the spin re- 
laxation rate, the spin is driven out of equilibrium. In the 
case of semi-integer S at zero applied magnetic field, this 
can result in the occupation of one of the two decohered 
ground states, g±^ and the depletion of the other, giving 
rise to a net magnetic moment {Sz)^ according to Eq. 
([5]). The population transfer takes place mainly through 
inelastic excitation of the spin from the ground state dou- 
blet Sz = ±5* to the first excited doublet, via spin- flip 
exchange. The transition rate where an ^ (majority) elec- 
tron from the tip spin-flips and goes to the surface reads 
(positive applied voltage in our sign convention) [18|] : 



7^^(|A + eF|)|(5+|5+|x+)p, 



(9) 



where we have assumed that \eV\ ^ A, ksT while 
refers to the excited state connected to In fact, the 
efficiency of the process is greatly enhanced when either 
bias or temperature are higher than the inelastic excita- 
tion energy, A 2^ (25* — 1)\D\ for half-integer spin S. 

In the case of integer spins, inelastic excitations also 
transfer population between the two tunnel-split ground 
states but, as the expectation value of the magnetic mo- 
ment in Eq. ([5]) at zero applied field in both states is 
null, {Sz) = 0, no matter which non-equilibrium distri- 
bution is achieved. In Fig. [Ih),d) we plot (5'^), defined 
in Eq. (j5j), as a function of a magnetic field for 3 situa- 
tions: zero bias, +10meV and -lOmeV, for both Fe and 
Mn on CU2N with finite tip polarization. At zero bias, 
we obtain the equilibrium Brillouin curve [19]. At finite 
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FIG. 1: (Color online) Schematic evolution of the energy spec- 
trum and degeneracies versus E for: a) integer spin S = 2 and 

b) half-integer spin S = 5/2. Magnetization curves of the Fe 

c) and Mn d) adatoms probed with a spin polarized tip with 
Vt = 0.74 for three applied bias: V = (thick black line), 
V = —10 meV (red line), and V = +10 meV (blue line). Here 
T = 0.5K, Tj/To - 0.5, Vr^ = 1, and pr = ps- 



bias spin transfer favors spin alignment parallel {V < 0) 
or antiparallel {V > 0) to the magnetic moment of the 
tip. The striking difference between integer and semi- 
integer spin is apparent in the figure. For integer spin, 
the magnetic moment is always null at zero field and the 
effect of bias is to heat the atomic spin decreasing the 
absolute value of (Sz) with respect to the zero bias case. 
For semi-integer spin, the atomic spin takes a bias de- 
pendent value at zero field. Thus, we find that current 
driven control of the magnetic moment of a single spin is 
only possible for semi-integer S. 

We now address the problem of back-action and the 
conditions under which a SP-STM can perform the quan- 
tized spin readout without perturbing the atomic spin 
state, avoiding the loss of the classical information. In 
other words, we look for a quantum non-demolition 
measurement [sl of the atomic spin using SP-STM, with 
the caveat that the atomic spin is decohered. The mag- 
netoresistive read-out [Eq. (j4])] is made possible by the 
tunneling exchange coupling between the quantum spin 
and the transport electrons. Specifically, it is based on 
the non-spin flip or Ising coupling, SzCTz^ which does not 
flip the atomic spin. However, if tunneling exchange is 
spin-rotational invariant, Eq. (j2j), the Ising term goes to- 
gether with the flip-flop terms, S~^(J~ + S~(J~^^ which 
induces atomic spin scattering with the selection rule 
ASz = ±1 and are responsible of the recording (spin- 
transfer torque). Thus, as in many other instances, the 
reading mechanism entails some degree of back-action on 
the probed system. The back-action occurs via inelastic 
spin- flip events, whose rate Finei takes off when either 
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FIG. 2: (Color online) a) Substrate mediated elastic spin re- 
laxation rate Ff^"^ in units of 70 = 7^'^(lmeV) for an ideal 
half-integer spin systems with D = — = — ImeV and 
T = 0.4K. b) Total relaxation rate {Finei + Tela) for S = 5/2 
at three different temperatures. Here Vt — —1, vs — 1, 
VT = 0.7, Tj/Tq = 0.5 and p|Tj = 0.01. 



mediated processes, Eqs. (j6]). 

In summary, we have studied the limitiations imposed 
by quantum mechanics to the use of quantum spins to 
store classical bits of information. We have found that 
classical information can be stored in quantum semi- 
integer spins, for which quantum spin tunneling is sup- 
pressed and quantum spin torque is possible, with uniax- 
ial anisotropy D <{). The storage time is limited, when 
un-observed, by the elastic spin- flip rate (Eq. [6|). Mag- 
netoresistive readout induces additional spin scattering 
given by the rate ([9]). Shot noise imposes a lower limit 
to the measuring time. Increasing using for instance 
few atom ferromagnetically coupled spin clusters, rises 
dramatically both the elastic and back-action lifetimes, 
as well as the decoherence rate facilitating the magnetic 
recording on a quantum spin. 
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and CONSOLIDER CSD2007-00010) and Generalitat 
Valenciana (ACOMP/2010/070). We acknowledge useful 
conversations with R. Aguado, G. Saenz and C. Untiedt. 



bias or temperature provides the excitation energy, and 
the elastic spin-tunneling assisted spin-flip, whose rate 
Tel depends only on ksT^ see Fig. [2j 

The condition for non-demolition readout is that the 
measuring time r is significantly shorter than the spin- 
lifetime, ^ (Einei + Eel). Rcgardlcss of the in- 
strumentation, the measuring time has a fundamental 
limit given by the condition that shot noise 51 should be 
smaller than the current contrast. A/ = AGV. For Pois- 

sonian noise, we have 51 = where / is the average 

current measured during r. If we define the average time 
for a single electron passage , Tg = e//, then, the limit 
imposed by shot noise is r ^ Tg. In other words, many 
tunneling events are necessary to perform the magneto- 
resistive single spin readout. 

Current experiments are done with / in the range 
of nA, which yields Tg ~ 0.2 ns, so that the measur- 
ing time is bound by below, due to shot noise, by Ins. 
State of the art instrumentation requires much larger 
measuring times. For instance, the use of lock-in in- 
troduces a more stringent bound to r, in the range of 
1/is-lms [IS EH. In Fig. Ill) we plot the substrate me- 



diated elastic spin relaxation rate Ff^*^ for and ideal spin 
system with D -^\E\ -ImeV and T 0.4K. 
This relaxation time grows exponentially with the spin 
S. For an experimentally sensible zero bias conductance 
G(0) ~ 0.01Go[2] {psTj = 0.1), relaxation time above 
1/is can be found in system with S > 5/2 with excitation 
energies A > 4 meV at T = 0.4i^. The bias dependence 
of the total relaxation rate is shown in Fig. [2]3)[18,]. In 
order to realize a non-demolition measurement, the bias 
should be kept |eF| <C A. As shown in Fig. [2]3), relax- 
ation rate in this regime is dominated by the substrate 
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